Anisotropic Dirac electronic structures of ^MnBi 2 (A=Ca and Sr) 



Geunsik Lee, 1 Muhammad A. Farhan, 1 Jim Sung Kim, 2 and Ji Hoon Shim 1, 2 ' 3, [j] 

'Department of Chemistry, Pohang University of Science and Technology, 
San 31, Hyojadong, Namgu, Pohang 790-784, Republic of Korea 
2 Department of Physics, Pohang University of Science and Technology, 
San 31, Hyojadong, Namgu, Pohang 790-784, Republic of Korea 
'' Devision of Advanced Nuclear Engineering, Pohang University of Science and Technology, 
San 31, Hyojadong, Namgu, Pohang 790-784, Republic of Korea 
(Dated: January 8, 2013) 

Low energy electronic structures in ^4MnBi2 (A=alkaline earths) are investigated using a first- 
principles calculation and a tight-binding method. An anisotropic Dirac dispersion is induced by the 
checkerboard arrangement of alkaline earth atoms above and below the Bi square net in j4MnBi2. 
SrMnBi2 and CaMnBi2 have different nature of the Dirac dispersion due to different stacking of 
nearby A layers, where the eclipsed- and staggered-type stacking cause two-fold and four-fold ro- 
tational symmetries, respectively. Using the tight binding analysis, we show the chirality of the 
anisotropic Dirac electrons as well as the sizable spin orbit coupling effect in a Bi square net. We 
suggest that Bi square net provides a platform for the interplay between anisotropic Dirac electrons 
and neighboring environment such as magnetism and structural changes. 

PACS numbers: 71.20.Ps, 73.90.+f 



I. INTRODUCTION 

Low energy electrons obeying the relativistic Dirac 
equation have been reported in the so-called Dirac mate- 
rials such as grapheneP^l Many intriguing physical prop- 
erties of Dirac fermions have been investigated both 
theoretically and experimentally, including unconven- 
tional quantum Hall effect, Klein tunneling, suppressed 
back scattering, and so onP Dirac materials are re- 
ported in topological insulatorj^Hsl j r0 n pnictidesp^ or- 
ganic conductor^ and inverse perovskitej^ and VO2- 
T1O2 interface pi and the list is rapidly growing these 
days. The effective Hamiltonian in all those systems is 
characterized by Pauli matrice, which is essential to man- 
ifest spinor-related phenomena. 

Recently an anisotropic Dirac cone in SrMnBi2 was 
found by density functional theory (DFT) simulation and 
confirmed by angle revolved photoemission spectroscopy 
and quantum oscillations P^l It is characterized by linear 
dispersions with strong anisotropy in the momentum- 
dependent Fermi velocity. Such Dirac band near the 
Fermi level is mainly from the Bi square net layerj^ 
and responsible for a finite Berry phase in Shubnikov- 
de Haas oscillations^ and the linear-field dependence of 
magnetoresi stanc eP^ Similar behaviors were also found 
in CaMnBi2) 15 * 16 ^This indicates that the low energy elec- 
tron motion in the Bi square net can be described by the 
Dirac Hamiltonian. So far, however, the origin of Dirac 
fermion in the Bi squre net has not yet been investigated 
clearly. So, this paper is devoted to study the origin of 
the anisotropic Dirac electron observed in SrMnBi 2 and 
CaMnBi 2 . 

This paper is organized as follows. Section II explains 
the crystal structure of two compounds and DFT calcu- 
lation details. In section III, the DFT band structures 
containing the Dirac dispersion are presented. In section 



IV, we clarify the mechanism of the anisotropic Dirac 
dispersion as well as the chirality of Dirac electron by 
using the tight binding (TB) method. Also the effect 
of spin orbit coupling (SOC) of heavy Bi atoms will be 
investigated. Finally we conclude our work in section V. 



II. SIMULATION METHOD 

Figure 1 shows crystal structures of CaMnBi2 (SG 129, 
P4/nmm) and SrMnBi 2 (SG 139, I4/mmm). Both com- 
pounds contain a square net of Bi atoms, which is indi- 
cated by the red (dark gray) balls. Ca or Sr atoms are 



P4/mmn 



W/mmin 



U AM 



Mii-Bi(2) 




Ca <d 

Bi(lfl 
Ca % 

(a) CnMnBi; 




Mn.-Bi(2> 



(b) SiMnBi; 



FIG. 1: (color online) Crystal structures of (a) CaMnBi2, 
(b) SrMnBi2, whose crystal symmetries and atomic positions 
are listed in Table I. Bi(l) or Bi(2) indicates the first or second 
type of Bi atoms. The orange (bright gray) line indicates the 
conventional unit cell. 
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located at the pyramid top of four base Bi atoms, where 
the vertical distance varys from 2.5 to 2.7 A. Each Ca 
or Sr layer has a checkerboard ordering with respect to 
the Bi square net. There exist another buffer layers con- 
taining Mn-Bi tetrahedrons. So there are two types of Bi 
atoms in the unit cell: Bi(l) in the square net and Bi(2) 
in the Mn-Bi layer. 

The stacking configuration of Ca or Sr layers above 
and below the Bi square net is different in these com- 
pounds. The neighboring Ca layers have the staggered 
lateral positions, while the Sr layers are eclipsed. Due 
to such difference, CaMnBi 2 (SrMnBi 2 ) has a primitive 
(body-centered) tetragonal Bravais lattice. Experimen- 
tal structural parameters of CaMnBi2 and SrMnBi2 are 
listed in Table I. 

Band structure calculations are performed with the 
full-potential linearized augmented plane-wave method 
implemented in the WIEK2K packaged The general- 
ized gradient approximation (GGA) by Perdew-Burke- 
Ernzcrhof (PBE) is used for the exchange-correlation 
potential.^ The radii of muffin-tin are set to 2.5 a.u. 
for all atoms. The number of k points used in the full 
Brillouin zone is 1000. Due to the magnetic Mn ions, the 
spin polarized calculation is carried out. Also the SOC 
is considered due to the presence of heavy Bi atom. 



III. DFT RESULTS 

In j4MnBi 2 (^=Ca and Sr), Mn 2+ has a 3d 5 electron 
configuration. The calculated spin magnetic moment 
within muffin-tin sphere is about 4 fig for both com- 
pounds, which show negligible changes under different 
magnetic configurations of Mn spins. In order to find the 
magnetic ground state, the total energies are calculated 
for various magnetic configuration such as ferromagnetic 
(FM), checkerboard- type antiferromagnetic (cAFM), and 
the stripe-type antiferromagnetic (sAFM) configurations. 
Their relative energies are listed in Table II. The ground 
state is found to be cAFM for both CaMnBi2 and 
SrMnBi2. The relative stability listed in Table II is qual- 
itatively consistent with the literature.^ The associated 
cAFM transition will be related to the magnetic transi- 
tion observed at ~ 290 K in the experiemntsl-^ Under 
the cAFM order, the interlayer exchange interaction is 
so weak due to large MnBi interlayer distance, and its 
contribution is as small as the order of magnitude of the 
energy tolerance (~ 1 meV/f.u.). In our calculation, we 
assumed the antiferromagnetic interlayer ordering of the 
Mn spins along the z axis, which does not alter the main 
feature of the band structures. The total energy differ- 
ence in Table II remains almost the same with considering 
the SOC. 

As shown in Fig. 1, two Sr (or Mn-Bi(2)) layers ad- 
jacent to the Bi(I) square net in SrMnBi2 are eclipsed 
when it is projected along the z axis, while the adjacent 
Ca layers are staggered in CaMnBi 2 . The total energy 
calculation results are consistent with the experimen- 
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tal observation that SrMnBi2 favors the eclipsed type, 
while CaMnBi2 favors the staggered type. So it is likely 
that the stacking type is mainly determined by the ionic 
size. Interestingly, another structurally related com- 
pound SrMnSb2 with the same stacking type as CaMnBi2 
is known to exhibit the chain-type reconstruction of Sb 
square net.^ However, CaMnBi2 does not show any dis- 
tortion in the Bi square net. The absence of such distor- 
tion in the Bi square net in CaMnBi2 is mainly due to 
the SOC. In fact we found that soft phonon modes of 42 
cm -1 (5.3 meV) lead to a chain-type reconstruction in 
CaMnBi2 without considering the SOC, but they disap- 
pear when the SOC is included. Thus strong SOC has 
an important role in the suppression of the chain-type 
distortion in CaMnBi 2 . 

Figure 2 shows the calculated band structures of 
SrMnBi 2 and CaMnBi 2 either without or with the SOC. 
Both compounds show spin-polarized Mn-driven bands 
near —3.0 eV and 1.0 eV with majority and minority 
spins, respectively. As indicated by the size of points, the 
bands near the Fermi level arise mostly from the Bi(l) 
p orbitals which are weakly hybridized with Sr or Ca d 
orbitals. Near the Fermi level, two linear bands cross at 
a certain wavevector, k , along T to M. When the SOC is 
included, two linearly crossing bands change to a quasi- 
linear bands with a small SOC-induced gap of 0.05 eV. 
The SOC-induced gap is much larger at X point, which 
removes most of states near the X point from the Fermi 
level. This is because the electronic states have smaller 
difference in energy near the X point than along the T-M 
symmetry line, and thus the SOC splitting is more sen- 
sitive to the nonvanishing SOC interaction. This will be 
shown more clearly in the section IV. C based on the TB 
analysis and the perturbation theory. 



TABLE II: Calculated total energies of ferromagnetic (FM) 
and the stripe-type antiferromagnetic (sAFM) configura- 
tions of in-plane Mn spins with respect to the ground state 
checkerboard-type antiferromagnetic (cAFM) configuration. 
The unit is eV per formula unit. 





cAFM 


sAFM 


FM 


SrMnBi 2 


0.0 


0.09 


0.55 


CaMnBi 2 


0.0 


0.08 


0.29 



3 



(a) CaMnBi 2 
cAFM cAFM + SOC 



(a)SrMnBi, 




r~~z 



(b) SrMnBi z 
cAFM cAFM + SOC 




FIG. 2: (color online) Band structures of (a) CaMnBi2, (b) 
SrMnBi2. The checkerboard- type antiferromagnetic (cAFM) 
order is assumed, and the cAFM+SOC means the inclusion 
of spin orbit coupling (SOC). The Bi-p atomic orbital contri- 
bution is indicated by the size of points in the cAFM results. 



The linearly crossing bands in CaMnBi2 and SrMnBi2 
have been a scribe d to the observed Dirac fermion-like 
dispersionP^HlfiM] j n pjg 3^ we snow the energy dis- 
persion near the Dirac point at k without the SOC. 
Here k\ and k% are defined as k\ — (k x + k y )/\/2 and 
k 2 = {—k x + ky)/V%, where (k x , k y ) corresponds to (0, 0) 
and (it /a, it /a) at T and M points, respectively. As shown 
in Fig. 3(a), the anisotropy, i.e. the ratio of Fermi ve- 
locities along fci and fc 2 directions, is as high as upto 
50. Also the hole and electron bands touch at the Dirac 
point. In case of CaMnBi2, shown in Fig. 3(b), over- 
all features in the energy dispersion near the k point is 
similar to that of SrMnBi 2 . However the gap introduced 
by the hybridization with the states from alkaline earth 
atoms is rather different in CaMnBi2. For SrMnBi2, the 
zero-energy gap is found only at a point k between T 
and M, while it is found along the continuous line in the 




FIG. 3: (color online) Band energy surface (left hand side) 
and dispersion curve (right hand side) near the Dirac point 
for (a) SrMnBi2 , (b) CaMnBi2 . fci is a coordinate along T-M 
direction, and its perpendicular axis is fe, where they are 
defined as ki = (k x + k y )/y/2 and = (— k x + k y )/\/2. 
Note the displayed range of ki and &2.' 0.283 < ki < 0.297 
and —0.10 < k-z < 0.10 in unit of 2n/a with their ratio 
Afc2/Afci ~ 14. The cAFM configuration is assumed, where 
the cAFM+SOC result is also shown by the dashed line in 
the right hand side. 



momentum space for CaMnBi2. These reflect the impor- 
tant role of the arrangement of alkaline earth atoms with 
respect to the Bi square net, which will be discussed in 
detail in section IV. A. In addition, when the SOC is in- 
cluded, a small gap is introduced at the Dirac point as 
mentioned above. However, the essential feature of the 
anisotropic Dirac dispersion is maintained. 

Figure 4 shows the calculated Fermi surfaces by con- 
sidering the SOC. Near T, hole pockets exist in both 
compounds. Around X or along T-M, there are strongly 
anisotropic pockects caused by the Dirac bands for 
CaMnBi2 and SrMnBi2. Such Dirac fermions are sup- 
posed to dominate the transport property mainly be- 
cause of their high Fermi velocities. From our calcu- 
lation, the carrier type is electron and hole near the 
X point and along T-M, respectively. Since the hole 
pocket in SrMnBi2 is absent in CaMnBi2, the hole car- 
rier density is greater in SrMnBi2 than CaMnBi2. This 
may be related to an experimental result that the pos- 
itive (negative) thermopower is observed in SrMnBi2 
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FIG. 4: (color online) Calculated Fermi surfaces. The zone 
center is regarded as F point. 



(CaMnBi2). 22 Mea nwhile , whether the observed Dirac 
fermion in CaMnBi j 15 * 16 ! is related to the electron pocket 
near X is questionable, where we have seen that lowering 
the Fermi level produces the small hole pocket along T-M 
like SrMnBi 2 . 



IV. TB ANALYSIS 

Our ab initio band structures for SrMnBi2 and 
CaMnBi2 indicate that the Dirac- like dispersion is mainly 
caused by the Bi(l) p x , p y orbitals. Depending on the 
stacking of Sr or Ca, they show Dirac point or the contin- 
uous band crossing line. In order to understand the main 
mechanism for different band structures, we carried out 
the TB analysis on both CaMnBi 2 and SrMnBi 2 . We 
construct the TB hamiltonian of Bi square net with and 
without interaction with Sr or Ca atoms, and discuss the 
mechanism for the Dirac-like electronic structures. Also 
the chiral nature of the Dirac electron as well as the effect 
of the SOC are investigated. 



A. Bi square net with unit-cell doubling 

Figure 5 represents a Bi square lattice interacting with 
Sr or Ca atoms. Because of the unit cell doubling due 
to Sr or Ca atoms, there are two Bi atoms at (0,0,0) and 
(a/2, a/2,0) in the primitive unit cell with the lattice con- 
stant a. They are denoted as B\ A and Bi B , respectively. 
First, we study the pristine Bi square net without consid- 
ering Sr or Ca atoms. Using relevant p x and p y orbitals 
of each Bi atom as basis, one can construct 4x4 TB 
hamiltonian after diagonalizing in the k space. Within 
the nearest neighbor approximation, the hamiltonian has 
a following form. 
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Here e p is the onsite energy of p x and p y orbitals. We de- 
note the basis for each column vector as 



FIG. 5: (color online) Bi square net with unit- 

cell doubling by (a) Sr, (b) Ca. The red or 
dark gray (bright gray) balls mean Bi (Sr or Ca) 
atoms. Atomic coordinates in unit of lattice con- 
stant a are Bi A (0,0,0), Bi s (l/2, 1/2, 0), Sr A (0, 1/2, 1/2), 
Sr s (0,l/2,-l/2), Ca A (l/2, 0,1/2), Ca B (0, 1/2, -1/2). 



where 4> x (r) and 4> y (r) are the atomic p x and p y orbitals 
of the Bi 1 atom with i = A, B. Following Slater-Koster 
parametrization, 23 the hopping terms V a p with a, ft = x 
or y are given below. 



V, 



2(£i7r + t\a) cos(k x a/2) cos(k y a/2) 
V xy = V yx = 2(i llr - t la ) sm(k x a/2) sm(k y a/2) 



Here ti a and are the nearest neighbor hopping pa- 
rameters for the a and it bonds. 

One can express Hq in a simpler form by using 2x2 
identity matrix /, and the Pauli matrix o~ x . 



e p I 

Vxx-f Vxy&x 



The orbitals in two Bi types interact only through the 
off-diagonal blocks V XX I + V xy a x whose eigenvalues are 
V xx +"f orb V xy with the relative phase j orb = ±1 between 
p x and p y orbitals. Thus the eigenvalues of H can be 
written as e p + -f atom {V xx + j orb V xy ), where -f atom = ±1 
denotes relative phase between two Bi types. 

We estimated the values of t\ a , to have a good 
agreement with a DFT band structure as shown in Fig. 
6(a). Bi p x and p v driven DFT bands are well repro- 
duced by using t la = 2.0 eV, t llt = -0.5 eV. The TB 
band structure with e p — is shown in Fig. 6(b). The 
additional bands from p z obital crossing the Fermi level 
are not considered in this analysis. 

The linear crossing at the Fermi level is found in Fig. 
6(b) due to folding of p x and p y bands from two Bi atoms. 
Such degeneracy point appears at every k point along 
the continuous line as we have seen in Fig. 3(b). The FS 
shown in Fig. 6(c) represents such degeneracy. So the 
unit-cell doubling in the Bi square net makes conduction 
and valence bands touch each other on the whole Fermi 
surface. 

We now consider Sr atoms located at a(0, 1/2, ±1/2) 
in the unit cell, as shown in Fig. 5(a). According to 
the DFT results of SrMnBi 2 , mainly Sr d orbitals par- 
ticipate in the hybridization with the Bi p band near the 
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FIG. 6: (color online) Band structures of single layer of Bi 
square net. (a) The band structure by DFT method. The 
contribution of Bi 6p x , 6p y atomic orbitals are indicated by 
the size of points, (b) The TB band structure of p x and 
p y with hopping parameters ti a =2.0 eV, tin = —0.5 eV. (c) 
Fermi surface plot by the TB method. 
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FIG. 7: (color online) Band structures of the SrBi lattice, 
(a) DFT band structure with Bi 6p m , 6p y orbital contribu- 
tions, (b) TB band structure, (c) E(k x , k y ) plot near the Dirac 
point, (d) anisotropic Dirac dispersion showing bands along 
F-M direction (ko+Sk, ko + Sk) and its perpendicular direction 
(fco + 5k, ko — 5k). 



Fermi level. Major contribution of Sr d orbitals comes 
from z 2 (^20 %) and xz/yz (~50 %) orbitals, but there 
are also non-negligible contributions from x 2 — y 2 and xy 
(~30 %). So we consider all five d obritals for the TB 
hamiltonian. 

The Hamiltonian matrix of the SrBi lattice, whose 
dimemsion is 14, is given below involving Hq. 

( H VsrBi \ 
nSrBi = I T/t T 

Here ta denotes the onsite energy of Sr d orbital, and / 
is an 10x10 identity matrix. VsrBi, the hopping term 
between Sr d and Bi p orbitals, is 4x10 matrix listed in 
Table II. 

As shown in Fig. 7, the TB band structure (Fig. 7(b)) 
is compared to the DFT result (Fig. 7(a)). By choosing 
u a = 1.5 eV, u n — —0.5 eV, a qualitative agreement 
is obtained between two methods. Around the linear 
crossing point, denoted by the circle in Fig. 7(b), the 
energy surfaces of two bands are plotted in Fig. 7(c). One 
can see that the degeneracy along the band crossing line 
is lifted except the Dirac point, inducing an anisotropic 
Dirac cone as clearly shown in Fig. 7(c). The anisotropy 
of the momentum-dependent Fermi velocity is estimated 
to an order of 10, as shown in Fig. 7(d). 



In the description of CaMnBi2, Ca atoms are located at 
a(l/2, 0,1/2) and a(0, 1/2,-1/2), as shown in Fig. 5(b). 
Similarly, Ca d orbitals are hybridized with Bi p orbitals. 
So 14 x 14 TB hamiltonian HcaBi is given below with the 
onsite energy td of Ca d orbital. 

„ (Hq V C aBi \ 

tlCaBi — T/t T 

The matrix elements of VcaBi are listed in Table III. 

In Figs. 8(a) and 8(b), we compare the band struc- 
tures obtained by the DFT and TB methods, respec- 
tively. Qualitative agreement is obtained by using the 
TB parameter u a = 1.2 eV, — —0.5 eV. In contrast 
to the SrBi layer, the degeneracy along the band cross- 
ing line is not lifted as shown in Fig. 8(c). Fig. 8(d) 
clearly shows the almost zero gap along the continuous 
line (fcrj + ^fcj ko — dk). 

The distinct electronic structures between the SrBi and 
the CaBi lattices arises from the different arrangement 
of alkaline earth atoms with respect to the Bi square 
net. This can be understood using the first-order de- 
generate states perturbation theory. Here, single Bi 
square net layer is regarded as an unperturbed system. 
At each k with < k x ,k y < 7T /a, two bands touch at 
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TABLE III: Matrix elements of VsrBi with coefficients c\ — (v / 3/v / 2)m (T , C2 = VZu^, c 3 = (\/2/4)u ff - (v^/v^K, c 4 = 
(v^/4)u CT + (I /v / 2)it^-, where k^ = k x a/2 and Ky = k y a/2. We have used u CT = 1.5, u n — —0.5 for the a, tt bonds, respectively, 
between Sr d and Bi p orbitals, where the resulting band structure is given in Fig. 7. 
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TABLE IV: Matrix elements of VcaBi with coefficients ci = (y/3/y/2)u a , C2 = V2u n , C3 = (v / 2/4)m (T — (V3/V2)u vr , C4 = 
(■s/6/4)u„ + (1 /v2)u x , where k^ = k x a/2 and K a = k y a/2. We have used u CT = 1.2, « w = —0.5 for the a, tt bonds, respectively, 
between Ca d and Bi p orbitals, where the resulting band structure is given in Fig. 8. 
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the Fermi wavevectors as shown in Fig. 6(c). Their 
symmetries are given by r y orb = +1 and <y atom = ±1. 
The unperturbed eigenvalues are e?(k) = V xx + V xy , 
e°(k) = — V xx — V yy . The only condition for the band 
crossing, i.e. e?(k) = (k), is V xx + V yy = 0. This re- 
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suits in tan(fca;a/2) tan(fc y a/2) = {tx a + ti„)/(ti a — ti v ). 
Assuming ti„ — to simplify our problem, we ob- 
tain cos(k x a/2 + k y a/2) =0. As a result, the k line 
of the band crossing satisfies k x + k y = n/a. The 
Dirac point ko = (ko,ko) has fco = n/2a. Also, from 
7° rb = +1, the unperturbed eigenstates can be writ- 
ten as a linear combination of local orbitals such as 

(^(r-R-r 

R and t = (a/2, a/2, 0). So they are expressed as 
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where ^° tom — 1 and j^ tom = — 1 give the bonding, 
anti-bonding states, respectively. 

Here we consider the perturbation by the additional 
potential of Sr atom layers. Each Sr atomic potential 
is assumed an isotropic potential v(r) to give the total 
potential V(r) = X^R l '( r ~ R ^ ( a /2)y) summed over 
whole lattice pointes, as shown in Fig. 9(a). In order 
to get the energy shift due to the perturbation, we need 



to evaluate the matrix elements Vu 
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FIG. 8: (color online) Band structures of the CaBi lattice: 
the scheme is same as that in Fig. 7, so refer the caption of 
Fig. 7. 



By asumming short-ranged function </>^_j_^(r), we consider 
the overlap integrals upto the second nearest neighbors, 
which are defined as follows. 
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J 2 = / dr^ +v (0)F(r)^ + „(ax) 



dv^ +y (0)V{v)^ +y (ay) 



V\i = Iq + Ii [1 + cos(fc :c a + fc. y a)] 

+ Jl [C0s(fca;0) + COs(fc a o)] 

+ (^2 + Ji) [cos(k x a) + cos(fc y a)] 

V22 = Jo — A [1 + cos(fe x a + kyd)] 
— Ji [cos(fca;a) + cos(fc a o)] 
+ (I 2 + J2) [cos(fc :c a) + cos(fc y a)] 

V12 = V 2 \ = 2(7 2 - J2) [cos(fc a a) - cos(k x a)] 

Along the band crossing line with respect to the Dirac 
point k , that is k = i(| + 5, § — S), it becomes V\\ — 
V22 = Iq with the only nontrivial element V12 = V%\- 
So diagonalization of V gives following energy eigenvalue 
shifts. 

Se = Io± 2(L 2 — J2) [cos(fc y a) — cos^a)] 
= 4(J 2 - J 2 )sin<5~4(/ 2 - J 2 )S 

It is linear in 5. The above derivation is confirmed with a 
numerical calculaton, whose result is shown in Fig. 9(b). 

Therefore, an anisotropic potential V(r) with I2 7^ J 2 
gives rise to the anisotropic Dirac band, which was the 
case of SrBi in Figs. 3(a) and 7. In case of the CaBi case, 
the condition I2 — J 2 = causes a vanishing energy shift, 
which also explains Figs. 3(b) and 8. 



Note I2 7^ J2 for the SrBi lattice, as illustrated in Fig. 
9(a), but I2 = J2 for the CaBi case. After straight- 
forward procedure, one can obtain following matrix ele- 



(a) 



ix Hj '] 





-1 1 2(xl(r> 
cos(A- t/)-cos(A- a) 



FIG. 9: (color online) (a) Illustration of local orbitals related 
to the Dirac bands in < k x ,k y < it /a, where <f> x + y (r) = 



c (r) + </> !/ (r))/v / 2with 7 



+ 1. The gray circles indicate 



a perturbing potential V(r) due to Sr atoms. I2 and J2 are 
the coupling terms of V between the second-neighbor orbitals 
(see the text), (b) Numerically calculated energy eigenvalues 
of the Dirac bands along the band crossing line, where the 
Sr potential shown in (a) is included. A linear variation with 
the offset from the Dirac point k x — k y = n/2a is confirmed, 
as derived in the text, as long as I2 7^ Ji- Otherwise the 
degeneracy is not lifted like CaBi case. 



B. chiral property 

In the anisotropic Dirac band of the SrBi lattice, orig- 
inal two bands in the folded zone are associated with 
symmetric and antisymmetric combinations of two Bi 
sublattice orbitals. In other words, the eigenstate is 



6^ + 7° 



"Vv ) + Of + l orb 4>y ) on one side of the Dirac 



point, while it is + 7 orfc <) - faf + 7° rb <) 
other side. One symmetry (<y a *°™ 1 = l) needs to be 
changed continously toward another symmetry ^ atom = 
— 1) around the Dirac point. This leads to a specific chiral 
nature around the Dirac point, which is clarified below. 



the 



Simplifying expressions as 



+ 7 



orbj^A 

y 



and 



0x +l° rb( t ) y = 02, the eigenstate can be written as a su- 
perposition of 4>i, cj>2, in other words V'k(r) = d\4>i{v) + 
^202(r) = Mi,^)- We find a unitary matrix that trans- 
forms </>i±</>2 into a spinor state with < a z >= ±1. Such 
transformation T is 



T = 



1 

71 



1 -1 
1 1 



For a given V'k(r), one obtains the expectation values 
of a a as (cr a ) = (di, d 2 | cr a T \d\, <i 2 } with a — x,y,z. 
For a = x,y, they are evaluated as (a x ) — d\d\ — 
^2^2) (cz) = —d\d2 — did^. Introducing a vector filed 
(n x ,n z ) = ((<j x ) , (a z )), we present {n Xl n z ) in Fig. 10 
for the hole and electron bands of SrBi near the Dirac 
point ko = (±feoj ko). The energy contours are shown as 



(a) 



4 



-0.2 -0 1 0.0 0.1 



02 k (2n/a) 



an expected form of the effective Hamiltonian H eS is 

H cti = kiv lz a z + k 2 v 2x a x , 
where the anisotropy ratio v z \jv x2 is an order of fO. 

C. spin orbit coupling 



(b) around (k k ) 




election 




■0016 

0.2908 0.2910 02912 02914 0.2916 0.290 8 0.2910 2912 2914 2916 
k (2«/n) k (2-Kla) 



(o) around (-k , k ) 



hole 





-0016 -0008 0000 0008 0016 -0016 -0.008 0.000 0008 0.016 
k, k ilnla) 



FIG. 10: (color online) (a) definition of fci and axis, and 
the positon of two Dirac points in the first and second quad- 
rants, (b) and (c) show (n x ,n z ) = ({cr x } , (<t z )) characters 
of the hole/electron states near the Dirac points (fco,fco) and 
{-k ,k ). 



a function of k±,k 2 for clearity as defined in Fig. 10(a). 
As shown in Figs. 10(b) and 10(c), the vector rotates 
2ir along the closed k loop, i.e. winding number 1. From 
Fig. 10(b), the arrow rotates clockwise as one follows the 
contour line in conter-clockwise in both hole and electron 
states at (ko, ko). We have checked that the chirality is 
the same for two Dirac points k and — k . In contrast, 
at (— fcrj, kg) the chirality is opposite to the case of (ko, 
fc ), as shown in Fig. 10(c). Accordingly, the opposite 
chirality is found at (fcrj, — ko). 

Around the Dirac point, one can construct the effective 
Hamiltonian in terms of a x and a z , as long as t he lin ear 
dispersion and the chiral nature are concerned! 24 * 25 ! At 
k = k + 5k, 

a— 0,x,y ,z 



The spin orbit interaction is typically described by a 
following potential. 



Hso(r) 



1 1 dU(r) 
i 2 c 2 r dr 



L S 



With replacing the radial integration by an effective con- 
stant Ago , we consider Hso = ^SoL ' S. 

The TB band structure of single layer of Bi square net 
is calculated with an additional term Hso- The matrix 



elements of L ■ S need to be found in the 



basis, 



instead of usual spherical harmonics representation, Y™ . 
By using the relationship 4> x — (Y^ + Y-^ )/s/2, 4> y — 
i{Yx — Y[~ )/v2, <f> z = Yj , we express the Hso hi terms 
of following basis for each of Bi" 4 and Bi s , 



j, a , b \ 










^*,t /' 


*vh ) ' 


/ > 


£a / » 





lA,B 



where f (I) indicates the spin. With this basis set, we 
obtain 6x6 Hso f° r each atom. 



/ i 1 \ 
4 i 
-1 -i 
0-1 -i 
i i 




H 



so 



^so 



\ 1 -i 



0/ 



Hso is added in the 12 x 12 H which additionally 
takes into account p z orbital and spin degree of freedom. 
The calculatd TB band structure is shown by the black 
solid line in Fig. 11(a), where Xso — 0.6, € x>y = 0.0, 
e z = —0.5. Also we have used t\ z cos(k x a/2) cos(k y a/2) 
with t\ z = 2ti n = —1 for the 7r interaction between 
(j)^ and (j)f . Qualitative agreement with the result of 
the first principle calculation, shown in Fig. 11(b), is 
excellent. In particular the SOC-induced splitting does 
not occur at every k point where the p x ^ y band crosses 
the p z band. In order to explain this, we show the 
latom = — 1(+1) bands whose local orbital forms are 
(j> x+y + latom<t> x+v and <f>f + y a tom<f>z> b y the dashed 
(gray solid) line in Fig. 11(a) without considersing 
Hso- Along r to M direction, each <fi z band crosses the 
(f> x + y bands twice. The leading contribution of the per- 
turbation Hso is the second order, which is given by 



H 



(2) 

so 



y e 



it^Hso^ 



7(ei C 



(o) ro) 



The 



main factor is that 



^m > \Hso\4'n'' > ' S ) is n °t vanishing only 



where v Q is the velocity parameters. From Fig. 10(b), when j a tom values are equal each other between tpm and 



9 



<a> (b) 




FIG. 11: (color online) Band structures of single layer of 
Bi square net with SOC by (a) TB method, (b) WIEN2k 
package. In (a) the sublattice symmetry j a tom = +1 or — 1 
is indicated by the gray solid or dashed line when SOC is not 
considered. 



(a) (b) 




FIG. 12: (color online) (a) TB band structure of the SrBi 
lattice with SOC. (b) Magnification of the region indicated by 
the circle in (a): the wavefunction symmetry near the Dirac 
point and the size of gap (~ 0.05 eV) are also indicated. 



tpn ■ Owing to this, H S q causes the degeneracy breaking 
only at the point where two bands with the same j a tom 
crosses. 

We try to reproduce the WIEN2k band structure of 
SrMnBi 2 with the SOC, shown in Fig. 2(c). Here we 
used slightly modified TB parameters of e X)V — -0.6, 
e z = —1.3, and Aso = 0.5, but maintaining all the other 
parameters same as listed in Table II. Also we do not 
consider the interaction between Bi p z and Sr d for conve- 
nience. The calculated SrBi TB band structure is shown 
in Fig. 12(a). Overall shape is consistent with the re- 
sult of the first principle calculation with disregarding 
Mn related bands. Also the small energy gap 0.05 eV, in- 
dicated in Fig. 12(b), agrees well with Fig. 3(a). As the 
symmetry property of 7 a t OTO is indicated in Fig. 12(b), 
the chiral property is not be altered by the SOC. 
V. CONCLUSION 

By using ab initio DFT method, we have shown that 
the degeneracy originated from the folded Bi p bands is 
lifted for SrMnBi 2 except for the Dirac point, which re- 
sults in the anisotropic Dirac band. But, in CaMnBi2 
the band crossing occurs along the continuous line in the 
momentum space. From the TB analysis, we conclude 
that the anisotropic potential created by Sr atoms is a 
main factor for the anisotropic Dirac band. Our study in- 
dicates that the nature of Dirac dispersion is sensitive to 
the effective nearest neighbor interaction in Bi square net. 
So, it will be interesting to investigate the Dirac nature 
of Bi square net with further changes such as structural 
distortion and magnetism. 
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